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Abstract 

We find the counterintuitive result that electrons move in opposite direction 
to the free electron E x B-drift when subject to a two-dimensional periodic 
potential. We show that this phenomenon arises from chaotic channeling tra- 
jectories and by a subtle mechanism leads to a negative value of the Hall 
resistivity for small magnetic fields. The effect is present also in experimen- 
tally recorded Hall curves in antidot arrays on semiconductor heteroj unctions 

but so far has remained unexplained. 
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We have long become used to the fact, that the Hall resistance of extended systems grows 
linearly with the magnetic field B apart from quantum phenomena like the quantum Hall 
effect. It thus came as a surprise when for small B in antidot superlattices the Hall resistance 
remained quenched and even took on negative values although the current is carried 
by electrons and not by holes. In spatially constricted systems like ballistic micro junctions 
of crossed wires similar phenomena had previously been reported and were attributed 
to ballistic rebound trajectories, where after colliding with a confining wall the electron 
is ejected into the opposite lead Since this explanation rests upon the finiteness 

and smallness of the junction, it cannot account for large and extended systems. Thus 
the negative Hall resistance in antidot superlattices (with as many as 10^ lattice cells) has 
remained mysterious and one may have wondered whether it is not an artifact or an edge 
effect. 

In the present paper we describe a general mechanism whereby the presence of a periodic 
superlattice potential causes a majority of electrons to move in a direction opposite to the 
free electron E x B-drift. It gives rise to quenching and negative Hall resistivity as a bulk 
property. In the regime of small magnetic fields there is a long-time tail associated with 
chaotic channeling trajectories, which tend to follow one of the channels of the 2D potential 
as they are trapped close to a tube of regular motion in phase space. The escape from these 
channels after crossing 3 to 1000 or more lattice cells is predominantly in opposite direction 
to the Lorentz force. The latter phenomenon bears a partial resemblance to the rebound 
trajectories responsible for quenching in ballistic microjunctions, where in a sense one deals 
with a single unit cell of an antidot lattice. However, it is the combination of the long- 
time tail and the asymmetric escape that gives rise to the negative Hall effect in periodic 
media, as one can see from two arguments. Rebound trajectories can have a larger curvature 
(radius ~ a) and would occur up to much higher magnetic fields {B/Bq ^ 0.5 in Fig. 
where quenching would be predicted but has disappeared in our extended system. This is 
the regime of magnetoresistance (Qxx) oscillations due to cyclotron-type orbits where 
the Hall resistivity (gxy) exhibits other features (Fig. H), which are also in good quantitative 
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agreement with experiment. Secondly, we find that the negative Hall resistivity is absent 
for infinitely steep antidots (Sinai billiard) where only a few channeling trajectories exist, 
whereas the rebound mechanism would be compatible with steep walls. Furthermore the 
quenching mechanism is not an edge but a bulk effect as it shows up even in the Hall 
resistivity and not only in the resistance. 

In micro junctions and lateral surface super lattices the dynamics of electron wave packets 
approaches the classical limit as the Fermi wavelength is smaller than the typical length scale 
of the lateral structure. We model the antidot superlattice with lattice constant a by the 
continuous potential 0] (Fig. |l] right inset) 

U{x,y) = Uo[cos{nx/a) cos{ny /a)]^, (1) 

with /3 = 4. The prefactor Uq of the potential is chosen such that the ratio of the dot 
diameter at the Fermi energy ep to the distance of adjacent dots is one-third, as in the 
experiments |l| . The classical approximation for the dynamics of an electron wave packet in 
the potential U{x,y) and a perpendicular magnetic field B || is then given by 4 coupled 
differential equations, which in scaled variables (x — > x/a, y y/a, U U/ep, t t/tQ, 
with to = {sF/rna^y/'^) read 

X = Vx Vx = 2V2{B/Bo)vy - dU/dx (2) 
y = Vy Vy = -2V2{B/Bo)v, - dU/dy. 

A chaotic trajectory for small magnetic field B/Bq = 0.1 is shown in Fig. 0, where Bq = 
2(2m£:^)^/^/ (ea) corresponds with a free cyclotron radius of a/2 of an electron with effective 
mass m. The trajectory follows one of the channels in the potential interrupted by local 
chaotic motion. Thus the free cyclotron motion in a magnetic field is destroyed by the 
periodic potential and the electron is forced into long episodes of linear motion. This is 
visualized in more detail in the left inset of Fig. ^ 

We investigate the motion in phase space (x, y, v^, Vy) for i? 7^ and = 0, by means of 
Poincare surfaces of section (x, Vx) at ?/(modl) = for various initial conditions. In Fig. |^ for 
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Vy > (left) and Vy < (right), this section displays a sea of chaotic motion and an island of 
regular motion which represents a tube of trajectories channeling forever in the y-direction. 



It contains periodic and quasiperiodic trajectories lying on cylindrical KAM-surfaces |T0 
around a central trajectory in phase space. They arise from nonlinear resonances between 
the degrees of freedom and are represented by closed loops in the Poincare surface of section. 
In the neighbourhood of such a tube of regular channeling trajectories, a chaotic trajectory 
can be trapped in an infinite hierarchy of permeable barriers, so-called cantori (not visible 
in Fig. ^). This appears to be the origin of the long channeling episodes of the chaotic 
trajectory in Fig. |l|. They are associated with a long-time tail as in similar systems which 
we previously studied |TT|JT^ . Here they will turn out to be essential for the negative Hall 
effect. We note that the tubes of regular motion in both directions {vy > and Vy < 0) as 
well as surrounding chaotic channeling trajectories still exist in a weak electric field E || e^,. 
This means that many trajectories move in opposite direction to the free E x B-drift and 
qualitatively explains the quenching of the Hall effect. For a more detailed account we now 
apply a linear response theory including regular as well as chaotic trajectories. 

We therefore divide phase space into regions of chaotic and regular motion corresponding 
to Fig. and write the conductivity cr as 



cr 



Prcr'' + pcO-", (3) 



where Pr and Pc are the portions of phase space volume occupied by regular and chaotic 
trajectories, respectively. Drude-type considerations for the straight regular trajectories in 
the channels lead to longitudinal conductivities a^^ = = prUre'^ /2m, where n is the 
electron density and r is the mean elastic scattering time, and to vanishing transversal 
conductivities a^y = ay^ = 0. Combining this with the conductivity of chaotic trajectories 
as derived in Ref. using hnear response theory yields 

(Jij =Pr<ylj + (1 -Pr)crO e'r Cij{t) dt, (4) 

where Cij{t) = {vi{t)vj{0)) is the correlation function of the unperturbed chaotic trajecto- 
ries (i.e. neither perturbed by an electric field nor by impurity scattering ) and the factor 



exp(— t/r) describes impurity scattering The numerical analysis involves the determi- 
nation of Pr as the volume of the outermost invariant KAM-surface enclosing the regions of 
regular motion divided by the total phase space volume. As these regular regions are small, 
however, the main contribution to the conductivities calculated according to Eq. stems 
from chaotic trajectories, whose correlation functions Cij{t) were determined from numerical 
simulations. The resulting Hall resistivity Qxy (Fig. ^ shows a quenching of the Hall effect 
and even a negative Hall resistivity as in the experiment The negative Hall resistance 
in antidot superlattices thus is a bulk effect due to chaotic electron trajectories. 

We now analyze the dynamical features leading to a negative Hall effect at small magnetic 
fields. The sign of the Hall resistivity g^y = o'yx/io'l^ + a^^) ( with a^y = —CTyx) is given 
by the sign of ay^, which according to Eq. is determined essentially by the correlation 
function Cyxit) of the chaotic trajectories. Figure^ shows the time dependence of Cyxit) for 
B / Bq = 0.05 where a negative Hall resistivity was found. Its initial increase still follows the 
free electron case but for later times there is a strong negative peak and a negative algebraic 
long time tail. This tail is responsible for the negative Hall effect as it yields a dominance 
of negative contributions in the integral over time in Eq. (^). Note that the mean elastic 
scattering time is r ^ 25 in the experiment [Q]. For larger magnetic fields {B/Bq > 0.12, 
see dotted line in Fig. ^) the tail disappeares and Cyx{t) exhibits an oscillatory decay where 
positive contributions to progressively outweigh the negative ones. The oscillatory decay 
for B/Bq > 0.12 is caused by chaotic trajectories that roughly follow cyclotron-type orbits 
around 1 or 4 antidots for a limited time and were shown to be responsible for pronounced 
peaks in the magnetoresistance QxxiB) (Ref. and dotted line in the inset of Fig. |]). Many 
rebound trajectories would have a similar curvature and would be expected in this 

magnetic field range {B / Bq ^ 0.5), yet they do not lead to a quenched or negative Hall effect. 
Instead one finds a prominent Hall plateau, which corresponds well with the experimental 
results p|,p^,p!5[1. In the experiment of Ref. [0, e.g. Bq ^ 0.55T. 



The long-time tail stems from chaotic channeling trajectories, which exhibit long episodes 
of straight motion (Fig. |l|) and dominate the dynamics below B/Bq = 0.12. The distribution 



ipii) of channeling lengths decays algebraically oc where u fluctuates around z/ = 3. 
Long-time tails seem to be typical for chaotic Hamiltonian systems and are believed to be 
related to cantorus hierarchies (see e.g. Ref. |T^|Tl|) as described in the context of Fig. 
but their origins are still not understood deflnitively. In the present case, after moving in 
a channel e.g. along the x-axis for a while (at least for 3 unit cells) the trajectories tend to 
leave the channel in the direction opposite to the Lorentz force — ev x B (which points to 
their left). For example at B/Bq = 0.05 a majority of about 65% of the trajectories show 
this anomaly (see Fig. § inset) and thereby yield the negative sign of the long-time tail. 
Intuitively the experimental situation can be imagined by an applied current transported 
mainly by channeled (chaotic) electrons moving in positive x-direction. Since they tend to 
leave the channel in negative y-direction they build up a negative Hall voltage. 

The role of the chaotic channeling trajectories for the long-time tail and the negative 
Hall effect can be described in more detail. For times t ^ a/vp and small magnetic flelds 
the decaying correlation function Cyx{t) is determined by long channeling episodes and can 
be expressed analytically in terms of the length distribution ■?/'(£) by 
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Cy^t) ~ ^ J^^^ m [P+W -P-W] di, (5) 

where f^;' is the Fermi velocity, avp is the mean absolute value of the contribution to the 
correlation function of individual channeling episodes, and {£) is the mean length of these 
episodes. p+{i) and p-{i) are the conditional probabilities that a trajectory trapped in a 
channeling episode of length i leaves the channel in positive or negative direction with respect 
to the Lorentz force. The asymptotic decay of Cyxif) obtained form Eq. (^) is indicated by 
the dashed line in Fig. ^. Here a and p+(£) — P-{t) were obtained approximately from 
geometrical considerations, while for the length distribution il){€) we used the numerical 
data. The long-time tail of Cyx{t) is relevant for the negative Hall resistivity up to times 
t ^ T, i.e. i ~ 30. The dashed line in Fig. | demonstrates that it is caused by the distribution 
ipi^) of chaotic channeling lengths. 

We emphasize that the results presented here do not depend sensitively on a particular 
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choice for the periodic potential. While for an infinitely steep antidot potential the effects are 
absent, realistic potentials with a small degree of softness produce a negative tail in Cyx{t) 
(we have checked this e.g. for /5 = 20 in Eq. (|1]). For the sample used in Ref. /3 = 4 seems 
realistic 0). The strength of the tail determines whether Eq. (|^) yields a negative Hall 
resistivity or only quenching. Furthermore irregularities of the potential lead to a reduction 
of the tail. (For a distorted potential with ±5% disorder in the antidot positions we still 
found a negative Hall resistivity.) This might explain why in some experimental situations 
quenching and in others a negative Hall effect was reported [0,0. 

We acknowledge motivating discussions with R. R. Gerhardts, K. v. Klitzing, and 
D. Weiss. This work was supported by the Deutsche Forschungsgemeinschaft. RK thanks 
Walter Kohn for the hospitality in Santa Barbara and for partial support by the NSF under 
DMR90-01502. 
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FIGURES 

FIG. 1. A chaotic trajectory in the antidot potential (right inset) most of the time shows 
straight paths along a channel even though a perpendicular magnetic field {B/Bq = 0.1) is applied. 
A magnification of such a channeled trajectory in positive and negative y-direction (left inset) shows 
an interplay of left turns due to the magnetic field and of right turns due to defiections from the 
antidot potential. 

FIG. 2. Poincare surfaces of section at y(modl) = for t;^ > (left) and Vy < (right) 
for a weak magnetic field B/Bq = 0.1. The closed curves in the center represent tubes of regular 
motion intersecting the plane perpendicularly in both directions. The quenching of the Hall effect 
is related to the survival of these tubes in both directions {vy > and Vy < 0) when a weak electric 
field E II e^; is applied. 

FIG. 3. Hall resistivity Qxy versus magnetic field B/Bq showing a negative Hall effect for 
B/Bq < 0.12 calculated from Eq. (^. Above B/Bq = 0.12 the inset shows a prominent plateau 
in Qxy for 0.4 < B/Bq < 0.9 (solid line) caused by chaotic trajectories revolving around 1 or 4 
antidots for a while (the dotted line shows Qxx after Ref. ^ for comparison). The small plateau in 
Qxy and the shoulder in Qxx at B/Bq ~ 0.25 are due to chaotic trajectories enclosing 2 antidots. 
(The Qxx and Qxy scales are not identical.) 

FIG. 4. The velocity correlation function Cyx{t) for B/Bq = 0.05 (solid line) exhibits a large 
negative peak and a negative algebraic tail. Both are caused by channeled chaotic trajectories 
(inset), which tend to leave the channel in a direction opposite to the Lorentz force — ev x B and 
thus give a negative contribution to the correlation. According to Eq. (Q) the area enclosed under 
Cyx{t) up to the elastic scattering time r = 25 Q is a measure for the Hall conductivity. — Above 
B/Bq = 0.12 other chaotic trajectories dominate, which revolve around 1 or 4 antidots and thereby 
produce oscillations in Cyx{t) (dotted line for B/Bq = 0.15). The dashed line shows how the length 
distribution 'ip{t) according to Eq. (|5|) gives rise to the tail for B/Bq = 0.05. 
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